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Abstract 


This work analyzes the water-gas-shift reactor design as component of the CO clean-up system of the ethanol processor for H) production 
applied to PEM fuel cells. The WGS reactor constitutes the element of greater volume of the processor motivating its optimization. A model-based 
reactor optimization for different reactor configurations permits to obtain both designs for reducing volumes and optimal operating conditions. 
The heterogeneous model used allows computing the optimal reactor length and diameter, and the optimal catalyst particle diameter. The model 
computes the constraints required for catalyst, such as maximum and minimum operation temperature. The volume is sensitive to the CO outlet 
concentration. According to the required CO conversion it is necessary more than one reactor unit for the case study analyzed. When considering 
the insulating material, there exists an optimal thickness that affects the final volume and the design variables. These results are useful for estimating 
the minimum and relative sizes that allows conventional reactor technology. 

Several reactors configurations are analyzed in order to state the limiting values of the main design variables. Specially, insulation conditions 


are studied in detail to access minimum total volumes. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


During the last decade, there have been important advances in 
fuel cells technology. Fuel cells are being developed for applica- 
tions to electrical energy generation and co-generation systems 
(coupled heat and power) in both stationary and mobile systems. 
Since fuel cells are high-efficiency energy converter devices and 
due to their low polluting emission levels, they become more 
and more attractive as a power generation alternative, specially 
in transportation industry [1,2]. 

Direct hydrogen fuel cell systems are clearly the preferred 
fuel cell operation mode. Nevertheless, the absence of a network 
for hydrogen distribution, and the risk and difficulties associ- 
ated with hydrogen storage and transport makes it difficult a 
widespread use of this alternative. On the other hand, a suc- 
cessful implementation of fuel cells in the short term can be 
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performed by fuel reforming. A promising route consists of the 
steam reforming of alcohols, mainly methanol and ethanol. The 
possibility of using ethanol reforming for producing hydrogen 
as combustible for fuel cells has generated an increased motiva- 
tion for investigating the reforming process of alcohols [3]. The 
ethanol shows advantages over fuels derived from fossil sources, 
since it is arenewable source, and has a neutral effect on the car- 
bon dioxide emissions, and can be obtained by fermentation 
[4]. 

In order to generate a hydrogen rich gas stream in the fuel 
processor, a feeding stream consisting of ethanol and water 
is converted into a stream with high hydrogen content. Fig. 1 
schematizes the global steam reforming process, which is a 
mature technology for synthesis gas, ammonia and hydrogen 
production from natural gas. 

Proton exchange membrane fuel cell (PEMFC) demands a 
CO-free hydrogen stream for operation. Since CO is adsorbed 
on the catalyst surface causing catalyst poisoning, it is neces- 
sary to reach CO concentrations less than 10 ppm for preventing 
irreversible damage and to facilitate the electrochemical reac- 
tion on the Pt electrode. The conditioning of the gas stream 
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Nomenclature 

ay external catalyst surface area per unitary reactor 
volume (cm7!) 

Ay cross reactor area (cm?) 

Ci component concentration along the reactor 
(mol cm~?) 

Cp; molar heat capacity of i component 
(kJ mol! °C7!) 

Cpr fluid mass heat capacity (kJ g7! °C7!) 

Csi component concentration inside the catalyst pore 
(mol cm~?) 

c$ i component concentration at the external catalyst 
surface (mol cm7?) 

Dett.i effective diffusivity in porous media (cm? s7!) 

D; reactor tube diameter (cm) 

Dp catalyst particle diameter (cm) 

Dins outer insulation diameter (cm) 

€ins insulating material thickness (cm) 

et reactor tube thickness (cm) 

F; molar flow of component i (mol s7!) 

G mass flow per unitary reactor area (g cm~? s7!) 

he energy transfer coefficient in the film 
(kJ cm~? s7! °C!) 

hy energy transfer coefficient (kJ cm~? s7! °C!) 

AH; heat of reaction (kJ mol!) 

kg mass transfer coefficient in the film (cm s7!) 

L reactor length (cm) 

Mi molecular weigth of i component (g mol~!) 

P pressure (atm) 

dev convective heat flux (kJ cm~? s7!) 

qrad radiative heat flux (kJ cm~? s7!) 

Q™ heat exchanged between the reactor and its sur- 
roundings per unitary bed volume (kJ cm~? s7!) 

rj reaction rate (mol g-cat™! s7!) 

Tamb environment temperature (°C) 

Te fluid temperature (°C) 

Tins surface insulating material temperature (°C) 

Ts temperature inside the catalyst pore (°C) 

Te temperature at the catalyst surface (°C) 

Tw reactor wall temperature (°C) 

Tc; critical temperature for the i component (K) 

Tr; reduced temperature for the i component 

V volume (cm?) 

Ve; critical temperature for the i component 
(cm? mol7!) 

Xco CO conversion 


Greek letters 


stoichiometric coefficient of component i in reac- 
tion j 

bed porosity 

insulation surface emissivity 

gas thermal conductivity (kJ em~! s7! °C7!) 


Aeff effective conductivity inside a particle 
iden ‘sc 4) 

uf gas viscosity (gcm! s7!) 

pf gas phase density (g cm7?) 

Pp catalyst density (g-cat cm7?) 

o Stefan—Boltzmann constant 
(5.670 x 107! kJ K~4 cm~? s7!) 


generated from the steam reforming of ethanol is partially per- 
formed by the water-gas-shift reaction (WGS). Afterwards, a 
final CO reduction is performed by the preferential oxidation 
reactor (CO-PrOx), where the oxidation of carbon monoxide 
is the desirable reaction. Nevertheless, undesirable combustion 
of valuable hydrogen also takes place. The reformer, WGS and 
CO-PrOx reactors are the nucleus of the steam reforming pro- 
cess, and contribute largely to the total volume and weight of 
the whole system. 

For fuel cell applications, a compact, efficient and reliable 
fuel processor is desirable. Process synthesis and design tasks 
are similar to other industrial reforming processes. However, 
the production capacity level required by fuel processors for 
vehicles or other similar devices is lower compared to industrial 
processes. Thus, process units of small size and specific designs 
are required. 

Since the WGS reaction rate is slower than the other reactions 
involved in the steam reforming process, and is limited at high 
temperatures by the thermodynamic equilibrium [5], the WGS 
reactor is the largest and heaviest process component. 

The aim of this work is focused on investigating the WGS 
reactor as a process component of a fuel processor for applica- 
tions in PEM-type fuel cells, and showing how the reactor design 
using mathematical programming techniques allows computing 
both reduced volumes and optimal operation conditions. Dif- 
ferent reactor configurations, reactor size and relative sizes of 
the reactor components (insulating material, reactor tube and 
reactive bed) can be evaluated; moreover, process performance 
bottlenecks and opportunities for optimization can be identified 
and analyzed. This knowledge can be used for process improve- 
ments such as lower unit weight and costs, and higher global 
efficiency. 


2. Water-gas-shift reaction 


The water-gas-shift reaction is widely applied in industry. It 
is used in the ammonia synthesis and hydrogen production pro- 
cesses from hydrocarbons reforming. Carbon monoxide reacts 
with water steam according to: 


CO+H20 + CO. +H2 AASog = —41.1 kJ mol! 


In industry, the main objective of the WGS reaction is to increase 
and to adjust the H2/CO molar ratio in the synthesis gas, and to 
remove CO from the off-gas. 

The reaction is moderately exothermic and its equilibrium 
constant decreases with temperature; high conversions are 
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C,H,OH 
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CO<1% CO < 10 ppm 


O, or Air 


Fig. 1. Scheme of the ethanol processor. 


favored at low temperatures. The equilibrium is virtually insen- 
sitive to the system pressure. Under adiabatic conditions, conver- 
sion in a single catalytic bed is thermodynamically limited. As 
the reaction progresses the reaction heat increases and the oper- 
ation temperature limits the conversion. The thermodynamic 
limitation of the reaction can be overcome using two or more 
catalytic beds with intermediate cooling, instead of a single bed. 

Although the WGS reaction is catalyzed by numerous mate- 
rials, it is industrially performed in two adiabatic stages using 
two different catalysts with intermediate cooling. The first stage 
(high temperature stage HTS) uses an iron-based catalyst (Fe/Cr) 
operating at 300-450 °C, and it converts the largest amount of 
CO. Typical CO levels in the outlet stream of a single reactor 
range between 2% and 4%. The second reactor (low tempera- 
ture stage LTS) operates at lower temperatures, between 180 °C 
and 230°C, using a copper—zinc catalyst supported on alumina 
Cu/Zn/Al. The LTS stage is favored thermodynamically, and 
can achieve a residual CO concentration around 0.1—0.3%. The 
catalyst composition is usually CuO/ZnO/A1O3, being very sen- 
sitive to poisoning by chlorides and sulfides, and to sintering. In 
addition, both catalysts are pyrophorics, i.e., they spontaneously 
release heat when are exposed to air after activation, increasing 
the temperature. 

In last years, it is noticed a renewed interest for this reac- 
tion, motivated by the promising application for coupling the 
steam reforming of alcohols with fuel cells for power and heat 
generation systems. 

Zalc and Löffler [5] mention that WGS reaction is limited by 
its intrinsic catalytic activity; therefore, WGS reactor will take up 
the largest volume in the fuel processor. For a methane processor, 
the authors estimated a volume of 226 cm? kW~! assuming an 
optimal temperature profile and for an outlet CO composition 
of 1% molar. 

Using commercial copper-based catalysts, Campbell [6] 
points out that the reaction is strongly controlled by pore dif- 
fusion at 200°C, and that copper-based catalysts are prone to 
lose activity due to sintering phenomenon. Thus, the operating 
temperature range must be limited to 150-250 °C. 

The efforts for improving the WGS reactor performance 
have been focused in a wide spectrum of subjects, including 
the development of more active catalysts [7—9], experimen- 
tal and theoretical studies to formulate kinetic expressions for 
novel catalysts [10-12], and theoretical studies related to heat 
and mass transfer phenomena in the catalyst particle structure 
[13]. Research is also concerned with applying and evaluating 
non-conventional reactors such as microreactors [14], mono- 
lith reactors [15] or membrane reactors to carry out the WGS 
reaction [16]. 


Giunta at al. [17] apply the one-dimensional heterogeneous 
model for simulating a WGS reactor using a commercial 
Cu/Zn/Ba/Al,03 catalyst. Based on simulation outputs, the 
authors recommend adiabatic operation of the reactor. 

In this paper a heterogeneous catalytic reactor model used for 
simulation and optimization is presented. Three different reac- 
tor structures are modeled as optimization problems and solved 
as study cases; specifically, it is addressed the optimal design 
and operation of: (a) a single adiabatic reactor, (b) two adia- 
batic reactors with intermediate cooling, and (c) reactor plus its 
insulation system. In addition, the effect of catalyst deactivation, 
changes on production scale and pressure effect over the design 
are analyzed. 


2.1. Mathematical model 


In this work, the reactor design is performed based on the 
one-dimensional heterogeneous model. This model offers higher 
accuracy for reactor design [18]. The mass, momentum and 
energy balances for the fluid phase and balances for catalyst 
particles are represented by Eqs. (1)-(7): 


Reactor model 
Fluid phase 
Mass balance 
dF; 


gg T Arkga(Ci — C3) (1) 


BC.:z=0 F; = F? 
Energy balance 
dTe 
GCP = hyay(Ts — Ti) — Q™ (2) 
BC.:z=0 F=7f 
Momentum balance (pressure drop along the catalytic bed) 


dP G? (1-— s) 
dz prDp ep 


L= 
1.75 +42 res E 8) 
e 


B.C.:z=0; P= P? 
(Spherical) catalyst particle 
Mass balance 


Nrx 
ld ds; 
az | Detit —= ) — pp) “air (Cs, Ts) = 0 4 
z ral eff iS dé ) Panis s» Ts) (4) 
dCs j 
B.C. : £ = 0; a0 
dé 
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S 
Cs,i = C3 i 


p dCs. (5) 


=) Ce 20h Das 


dé |p,/2 


Energy balance 


Nrx 
ld aT: 
(ene?) + pp) (-AH)j;rj(Cs,Ts)=0 (6) 
j=l 


£2 dé dé 

dTs 
B.C. :£=0; —2=0 

d 

S 
pa os dT. (7) 
~ 2) AdTS — TÀ = Aet — 
dg Dp/2 


Eqs. (5) and (7) represent the mass and energy balances for the 
fluid-particle interphase, formulated as boundary conditions. In 
Eq. (2), Q™ is the heat exchanged between the reaction system 
and its surroundings. In the adiabatic operation mode Q™ =0. 
When insulation is incorporated into the model, this term eval- 
uates the heat transferred to the environment. 

Correlations for the fluid heat capacity (Cpr), viscosity (ur) 
and thermal conductivity (àf) are those given in [19,20] (see 
Appendix A). Correlations for the mass and energy transfer coef- 
ficients in the film (kg and hg, respectively), as well as expressions 
to estimate the effective diffusivity in porous media (Desf) are 
those given in [21]. Gas ideal behavior is assumed to evaluate 
the mixture density of the gas phase (pf). The average molecular 
weight of the gaseous mixture and the molar flow entering the 
reactor were used to compute the mass flow per unitary area (G). 
Heat transfer coefficients are those given in [22]. 

The bed porosity (€p) is estimated by the expression given in 
[22]: 


_ (D/D = 2° 


1 
(D,/Dp)* 


£b = 0.38 + 0.073 (8) 


For spherical shape particles, the external catalyst surface area 
per unitary reactor volume (ay) is computed by ay = 6(1—ép)/Dp. 

Summarizing, molar flow of component i (F;), fluid temper- 
ature (Tf) and pressure along the reactor bed (P) are calculated 
from differential Eqs. (1)-(3). Component concentration (C Š. D 


and temperature (T$) at the external catalyst surface and inside 
the catalyst pore (Cs,;, Ts) are obtained from Eqs. (4)-(7). 

A fully developed flow is considered at the reactor intake. 
However, it should be taken into account that the flow into the 
fixed-bed reactor is generally achieved by means of a feed pipe 
and a distribution hood. These must therefore be constructed so 
that the fixed bed is uniformly traversed, and the gas residence 
time in each flow path of the fixed bed is the same. The effects 
of the entrance length and the reactor hood design to achieve the 
uniform gas flow to the catalyst bed are not considered by the 
model. 


2.2. Optimization model 


The optimization problem is formulated to obtain the optimal 
operating conditions and equipment size aiming at minimiz- 
ing the system volume. Following, the optimization problem is 
formulated, i.e. the objective function, decision variables and 
constraints are specified. 


Min (Viotal) 
T: Li, Di, Dp, eins 


Objective function 
Decision variables 


Interior point constraints 
Catalyst temperature 
External insulator temperature 


150°C < Ts < 250°C 
Tins < 60°C 
Final point constraints 


Admissible pressure drop 
CO molar fraction 


(P® — PSY/P? < 0.3 
yco < 0.003 


Design constraints 


Plug flow condition L,/Dp >30; D/Dp >10 


As mentioned, it is intended to minimize the total volume of 
the system. In WGS reactor, the total volume is determined by 
the volume occupied by the catalyst bed and the volume occupied 
by the insulating material. The optimization problem determines 
the optimal reactor length (Z;), reactor diameter (D+), catalyst 
particle diameter (Dp), and insulating material thickness (eins) 
that minimize the total system volume. For adiabatic operation 
case, the constraints corresponding to the insulating material 
temperature (Tins) and insulation thickness (eins) as decision 
variable are not included into the model. 

Due to the exothermic nature of the WGS reaction it is 
desirable to keep the reactor temperature within a certain oper- 
ating range. An upper temperature bound is set for avoiding 
catalyst sintering. At low temperature, the reaction rate dimin- 
ishes, the kinetic expression is not longer valid, and water 
condensation occurs. Then, it is necessary to impose a lower 
temperature bound. The reactor inlet temperature (To) is a 
decision variable, and its value results from the optimization 
problem. 

The resulting partial differential algebraic equations (PDAEs) 
are implemented and solved using gPROMS (general Process 
Modeling System) [23]. gPROMS is a general purpose model- 
ing, simulation and optimization system software. PDAEs are 
systems of equations that, when are discretized in spatial direc- 
tions (using any of the traditional numerical methods), result in 
differential algebraic equations (DAEs). DAEs are essentially 
sets of ordinary differential equations (ODEs) where some of the 
variables are constrained by algebraic relations. Then, the result- 
ing initial value DAE problem is integrated in the z direction 
present in the system. 

The algorithm used in gPROMS for solving the differential 
equation system is based on a backward differentiation formula 
(BDF) type method. The composition and temperature profiles 
established inside the catalyst particle, are obtained by dis- 
cretizing Eqs. (6) and (7) using orthogonal collocation on finite 
element method. Finally, the optimization algorithm used is the 
single-shooting method, which is also available in gPROMS 
environment. 
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2.3. Kinetic model 


The reactor is optimized considering the most probable con- 
ditions for the WGS reaction in a small-scale ethanol reformer 
using a commercial Cu/ZnO/AI103 catalyst from Sud-Chemie. 
The kinetic expression used is that provided by Choi and Stenger 
[11]: 


47400 Peo, Pa, 
rco = 82.2 exp C) (FeoPis0 = oe) (9) 
eq 


where keq = Pco, Pu, /PH,0Pco, and can be computed from 
expressions given in [11]. 


3. Results and discussion 


If a fuel processor for generating 1 kW of power is assumed, 
33.3 mol h`! of hydrogen is required [11]. This production rate 
can be reached if about 10 mol h~! of ethanol is processed in the 
steam reformer operating at 700°C with a water/ethanol molar 
ratio of 4.0. The molar flow rate values of the reformed gases 
obtained at these operating conditions are the ones reported in 
[24]. These values are assumed to be the input specifications 
to the WGS unit. A heat exchanger is needed for reducing the 
gas temperature to the operating range for favoring the WGS 
reaction. Thus, the WGS reactor’s inlet temperature is a decision 
variable that results from the optimization problem. The molar 
flow rates and compositions entering to the WGS reactor are 
listed in Table 1. 


3.1. Adiabatic reactor 


Firstly, the adiabatic operation of a single-stage reactor is con- 
sidered. Three design goals for the outlet CO concentration are 
specified: 1, 0.7 and 0.3% for each case. Although the CO-PrOx 
reactor can operate with gaseous mixtures with CO concentra- 
tion ranging between 0.2% and 2.0%, a design goal of 0.3% may 
assure good system efficiency. 

Since hydrogen is generated by the WGS reaction but is unde- 
sirably consumed during the oxidation of CO in the CO-PrOx 
reactor, there exists a trade-off between the ethanol processor 
efficiency and the volumes of both reactors. Consequently, an 
integrated analysis considering both aspects is necessary to find 
out the suitable CO conversion target in the WGS unit. 

The heterogeneous model allows computing the main vari- 
ables of the system geometry (e.g. optimal reactor length and 
diameter, and optimal catalyst particle diameter) as well as the 
main operation variables—e.g. optimal inlet reactor temperature 
restricted to the thermal profile imposed for avoiding material 


Table 1 
Inlet molar flow rates and compositions to the WGS reactor 


Molar flow rate (mol ho!) Molar fraction 


Fon, 4.7 YCH, 0.057 
Fu, 36.0 YH 0.435 
Fco 6.6 yco 0.080 
Fo 26.4 Ymo 0.319 
Fco, 9.1 YCO, 0.110 


Table 2 

Optimization results for a single-stage adiabatic WGS reactor 

CO output (%) 0.30 0.70 1.00 
L (cm) 36.50 10.00 10.90 
D, (cm) 7.20 7.39 6.39 
Dp (cm) 0.057 0.057 0.05* 
T? °C) 127.1 164.8 168.7 
Volume (cm?) 1460 432 339 


à The calculated value corresponds to the lower bound. 


deterioration. The optimization results obtained are included in 
Table 2. 

It is known that there exists a relationship between the cata- 
lyst particle diameter and the pressure drop along the catalytic 
bed. As the catalyst particle diameter decreases, the mass trans- 
fer resistance inside the porous matrix of the catalyst decreases, 
whereas the system pressure drop increases along the catalytic 
bed. For almost all analyzed cases, the calculated particle diam- 
eter Dp is 0.05cm, which is the lower bound value assigned 
to this decision variable in the optimization problem. From the 
point of view of kinetics and diffusions aspects, the use of small 
particles favors the reactor performance. Then, the diameter of 
the cross section of the reactor should be large enough to fulfill 
the pressure drop constraint imposed. For all analyzed cases, 
the pressure drop is about 0.1% with respect to the inlet pres- 
sure (1 atm). This value is negligible compared to the maximum 
pressure drop usually accepted and specified in the optimization 
model (30%). 

An increase of the reactor temperature is thermodynamically 
unfavorable for the water-gas-shift reaction. In adiabatic opera- 
tion mode, for achieving a CO composition of 0.7% in the outlet 
reactor stream, an increase of the reactor volume about 30% 
is needed with respect to reach 1.0% of CO in the off-gases; 
whereas a reaction volume increment around 330% is needed 
to reach a CO composition of 0.3%. If a higher CO conversion 
is required, a lower inlet temperature is needed in order to keep 
the temperature profile along the reactor below the maximum 
permitted value. 

In order to achieve the goal of 0.3% of CO composition 
using a single-stage unit operating in adiabatic mode, it was 
necessary to diminish the lower bound imposed to the operation 
temperature to 127°C. It means that a reactor design that ful- 
fills all process, design and operation constraints imposed, is not 
feasible as the problem was originally formulated. 

Fig. 2 shows the simulation results corresponding to the opti- 
mal reactor design for a CO composition of 0.3% in the off-gas 
stream. The temperature and CO molar fraction profiles along 
the reactor axis are depicted. It can be observed that the reac- 
tor volumes to achieve 1.0% and 0.7% of CO concentration in 
the off-gas are larger compared to the volumes computed by 
optimization when those concentrations are specified as design 
constraints. 


3.2. Two adiabatic reactors in-series 


An alternative arrangement to achieve a final CO composition 
of 0.3% in the off-gas stream, keeping the reactor temperature 


0% 
- - -10% 
---- 50% 


Activity Loss 
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Fig. 2. Simulation results for the design of a final CO concentration of 0.3%. 
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Fig. 3. Optimization results with two adiabatic catalytic beds in-series. 


in the desired operation range, is to consider two catalytic beds 
in series with an intermediate cooling unit (Fig. 3). 

To formulate this problem, a new decision variable has to 
be considered: the inlet temperature to the second catalytic bed. 
Fig. 3 shows the results obtained from the optimization problem, 
while Fig. 4 shows the simulated temperature and composition 
profiles along the whole system. 


Temperature [°C] 
CO Molar Fraction 


Second Stage 


9 9, 
0.7% CO 


0 5 10 15 20 25 30 35 
Reactor Axis [cm] 


Fig. 4. Temperature and CO molar fraction profiles along the catalytic beds. 
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Fig. 5. CO molar fraction profile for different catalyst deactivation levels. 


The total volume computing both catalytic beds (847 cm?) 


is smaller compared to the volume resulting from a single unit 
(1462 cm). However, it should be noted that the volume occu- 
pied by the heat exchanger is not evaluated. The heat exchanger 
design will depend on the cooling fluid used which may be 
available from the process streams involved in the global pro- 
cess favoring the system energy integration, and resulting in an 
increased process efficiency. 


In addition, this reactor configuration model allowed finding a 


reactor design that reaches the desired CO level (0.3%) without 
diminishing the temperature below 150°C, value fixed as the 
lower operating limit. 


3.3. Catalyst deactivation 


Catalyst deactivation was not considered in the previous reac- 


tor configuration designs. The catalyst used showed an activity 
loss of 10% for a 250 h operation period [11]. Fig. 5 plots simula- 
tion results for different catalyst deactivation levels considering 
the design obtained for the adiabatic reactor to achieve a CO 
molar fraction of 0.7%. It can be observed that when the activity 
loss increases the reactor is not being able to fulfill the desired 
conversion. However, in practice, the effect of the activity loss 
can be compensated by increasing the inlet reactor temperature. 
Table 3 shows optimization results when catalyst deactivation is 
taken into account at the reactor design stage. 


Table 3 

Optimal reactor design considering catalyst deactivation 

Activity loss (%) 0 10 50 
CO output (%) 0.70 0.70 0.70 
L (cm) 10.00 10.00 14.84 
D (cm) 7.39 7.79 8.46 
Dp (cm) 0.05" 0.05" 0.05* 
T? °C) 164.8 164.9 164.7 
Volume (cm?) 432 476 833 


è Computed value is at the lower bound value. 
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Fig. 6. Scheme of the reactor+insulator system. 


3.4. Thermally insulated reactor 


In practice, adiabatic operation is performed by attaching an 
insulating material layer to the wall of the reactor tube for pre- 
venting heat loss to surroundings. To analyze its influence on the 
final size and reactor behavior, heat loss rate to the environment 
is included into the heterogeneous reactor model. 

For a tubular reactor where a reversible exothermic reaction 
takes place, it is known that the maximum conversion is obtained 
by a unique temperature profile that maximizes the reaction rate 
at every point along the reactor. Specifying an adequate extent 
of insulation represents a pragmatic approach to affect the opti- 
mal heat exchange, and thereby to approximate the optimum 
temperature profile [25]. 

The insulating material used in the model consists of a layer 
of calcium silicate covered with a thin layer of aluminum. The 
heat loss rate to surroundings is modeled as an energy trans- 
fer phenomenon by heat conduction through materials involved, 
and by convection and radiation from material surface to envi- 
ronment (see Appendix B). The insulating material conductivity 
is considered as a quadratic function of temperature according 
to the correlation k(T)=ax+b,y T+cx TŽ. Heat loss by con- 
vection and radiation are evaluated by expressions given by 
Koenig [26]. Three insulations can be localized; side insula- 
tion, which depends on reactor wall temperature; front and end 
insulations, which depend on input and output reactor tem- 
peratures, respectively (see Fig. 6). Front and end-insulation 
were calculated considering the whole surface determined 
by the reactor-insulation outer diameter. Consequently, pos- 
sible effects due to the feed pipe to the reactor were not 
considered. 


Owing to that the thickness of commercial-steels tube 
increases with the diameter, it is necessary to model this behavior 
to penalize this effect on the total volume. A quadratic expression 
of the thickness tube as function of the internal diameter is con- 
sidered (e; = a + bD; + cD?). This correlation was obtained by 
fitting data from commercial tubes. 

The optimization variable set considered involves the same 
variables as in the adiabatic case plus the thickness of the insulat- 
ing material to minimize the total volume of the system (reactor 
plus insulator), keeping a suitable outer wall temperature. 

Firstly, a case with only side insulation is considered. Table 4 
summarizes the optimization results for different outer wall tem- 
perature values. The insulator thickness is a key decision variable 
that affects the equipment size. If the outer wall temperature Tins 
is restricted to values lower than 60°C, the total volume of the 
system is 1760 cm. Although the total volume (catalytic bed 
plus insulator) is larger than that computed for a single adia- 
batic reactor, the volume of the reactive catalytic bed itself is 
smaller. 

Fig. 7 shows temperature and CO molar fraction profiles for 
two design cases (Tins < 50°C and Tins < 120°C). The total heat 
loss in each case is about 15 W and 60 W, respectively. If a higher 
insulator wall surface temperature is allowed, the total volume 
is reduced to 955 cm°. It can be observed from Table 4 that 
for the cases where the outer temperature is limited to values 
higher than 80°C there is no need for an insulating material, 
i.e. the calculated insulator thickness value is at its lower bound 
(0.001 cm) and the outer reactor metallic tube wall reaches that 
temperature. 

These results indicate, as it is known, that losing heat to sur- 
roundings favors the equipment design, i.e. a reduced size is 


Table 4 

Optimization results for the reactor—insulator system (side insulation only) 

Tins (°C) <50.00 <60.00 <80.00 <120.00 
CO Output (%) 0.30 0.30 0.30 0.30 

Ly (cm) 18.60 13,31 15.56 25.36 

D; (cm) 9.66 11.00 9.00 6.00 
Dp (cm) 0.05* 0.05* 0.05* 0.05 
€ins (cm) 1.72 0.37 0.001 0.001? 
T? (°C) 148.5 148.7 154.9 187.5 
Bed vol. (cm?) 1360 (46%) 1260 (72%) 990 (79%) 720 (75%) 
Tube vol. (cm?) 350 (12%) 310(17%) 270 (21%) 240 (25%) 
Ins. vol. (em?) 1260 (42%) 190 (11%) 0(0%) 0(0%) 
Total vol. (em?) 2970 1760 1260 960 


* Computed value is at the lower bound value. 
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Fig. 7. Temperature and CO molar fraction reactor profiles: (a) Tins < 50°C and 
(b) Tins < 120°C. 


computed because the reaction is favored by thermodynamic 
effects. Moreover, a small heat loss diminishes the reactive bed 
volume compared to an adiabatic reactor. 

The effect of considering front and end insulation over 
the system design is shown in Table 5. The thickness of 
these insulation sections depend on the inlet and outlet reac- 
tor temperature. When the problem considers these sections 
the optimization results vary slightly, showing smaller reactor 
diameters. 

In order to verify the validity of the 1-D heterogeneous 
model used at the design stage (optimization model), simula- 
tion results are compared to simulation outputs obtained using 
a 2-D pseudo homogenous model. Fig. 8 depicts the tem- 
perature and CO molar composition profiles for the reactor 
design obtained adopting a maximum outer wall temperature 
of 50°C. The difference between the radial average values 
from the 2-D model and those predicted by the 1-D model is 
negligible. This agrees with the analysis performed by Koven- 
klioglu and DeLancey applied to a SO2 converter [25]. By 
means of a parametric analysis of the Nusselt number and 
using simplified models they showed that both 1-D and 2- 
D models predict approximately the same optimum thickness 
value. 
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Fig. 8. Comparison between 1-D heterogeneous model and 2-D pseudo homoge- 
nous model (radial average values). 


3.5. Effect of the change on the production scale 


Table 6 lists the optimization results for larger production 
scale of the processor to reach different power generation targets. 

It is worth to note that the model predicts an increase in 
the particle diameter when increasing the hydrogen production 
(25 kW and 50 kW cases). As the flow rate of the reactive mixture 
increases, a particle size increment is computed to maintain an 
admissible pressure drop through the reactor. 

These results indicate that a single reactor processing a flow 
rate for a power generation target of 50 kW would need for a vol- 
ume 70% larger than five reactors in parallel of 10 kW each. This 
behavior is owed to the diffusional resistance when increasing 
the particle size, decreasing the catalyst effectiveness as can be 
seen in Fig. 9. The effectiveness factor along the reactor axis for 
different designs is depicted in Fig. 9. Table 6 shows the activity 
loss due to diffusional resistance, the catalyst volume (bed vol- 
ume) per unit of kW generated increases with the target power. 


Power Target 


Effectiveness Factor 


0.0 0.2 0.4 0.6 0.8 1.0 
z/Lt [Dimensionless] 


Fig. 9. Effectiveness factor vs. dimensionless reactor length. 
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Table 5 
Optimization results for the reactor—insulator system (side, front and end insulations) 
Tins (°C) <50.00 <60.00 <80.00 <120.00 
CO output (%) 0.30 0.30 0.30 0.30 
Ly (cm) 22.23 15.02 15.71 25.54 
D; (cm) 8.14 9.82 8.95 5.99 
Dp (cm) 0.05" 0.05* 0.05* 0.05* 
eins (cm) Side 2.02 0.63 0.001? 0.001° 
eins (cm) Front 2.22 1.31 0.64 0.32 
eins (cm) End 4.68 2.88 1.34 0.40 
T (°C) 148.7 148.5 154.3 187.2 
Bed vol. (cm?) 1160 (29%) 1140 (50%) 990 (70%) 720 (73%) 
Tube vol. (cm?) 330 (08%) 290 (13%) 270 (19%) 240 (24%) 
Ins. vol. (cm?) 2530 (63%) 840 (37%) 160 (11%) 30 (03%) 
Total vol. (em?) 4020 2270 1420 990 

* Computed value is at the lower bound value. 

Table 6 
Optimal reactor designs for different power generation targets (side insulation only) 
Power target (kW) 1 10 25 50 
Tins CC) <50.00 <50.00 <50.00 <50.00 
CO output (%molar) 0.30 0.30 0.30 0.30 
Ly (cm) 18.60 105.58 207.23 317.75 
D; (cm) 9.66 13.40 17.75 24.28 
Dp (cm) 0.057 0.05* 0.11 0.16 
ins (cm) 1.72 1.65 2.10 2:32 
T? °C) 148.5 149.1 149.2 149.2 
Bed vol. (cm?) 1260 (72%) 14880 (55%) 51260 (57%) 147180 (63%) 
Tube vol. (cm?) 310(17%) 3240 (12%) 9270 (10%) 19990 (09%) 
Ins. vol (cm?) 190 (11%) 8990 (33%) 29190 (33%) 65410 (28%) 
Total vol. (em?) 1770 27100 89710 232570 
Bed vol. per kW 1260 1490 2050 2940 
Total vol. per kW 1770 2710 3590 4650 


à Computed value is at the lower bound value. 


In conclusion, it is convenient to combine reactors in parallel 
using smaller particle diameters. 

For two reactors with intercooler the results do not show a 
particle increment (Table 7). The temperature decrease between 
reactors allows achieving a design that satisfies the admissible 
pressure drop and temperature constraints without increasing 
the particle diameter. 


3.6. Pressure effect on reactor design 


Previous results were obtained considering | atm as pres- 
sure design because of kinetic expression was determined 
for that condition [12]. As commented by Giunta et al. [16] 
the CO conversion is favored when pressure increases. The 
kinetic expression used in this work depends quadratically 
with pressure. An operating pressure increment favors the 


Table 7 

Two-stage reactor design with insulations for 1 kW and 50 kW 

Power target (kW) 1 50 

Tins (°C) <50.00 <50.00 

CO output (%) 0.30 0.30 

Reactor stage #1 #2 #1 #2 

L (em) 5.97 6.24 20.07 25.13 

D; (cm) 13.88 6.36 39.00 29.31 

Dp (cm) 0.057 0.057 0.05 0.057 

eins (cm) Side 0.58 1.24 2.26 2.14 

€ins (cm) Front 2.27. 2.58 3.25 4.37 

eins (cm) End 5.35 2.46 7.18 4.47 

T? (°C) 148.00 164.80 164.63 195.74 

Bed vol. (cm?) 900 (31%) 200 (23%) 23970 (47%) 16950 (50%) 
Tube vol. (cm?) 190 (07%) 60 (07%) 3150 (06%) 2390 (07%) 
Ins. vol. (cm?) 1790 (62%) 590 (69%) 68110 (46%) 14420 (43%) 
Total vol. (cm?) 2880 (77%) 850 (23%) 47330 (60%) 33760 (40%) 
Total vol. #1+ #2 3730 cm? 84510 cm? 


à Computed value is at the lower bound value. 


476 J.A. Francesconi et al. / Journal of Power Sources 173 (2007) 467-477 


reaction rate and the required bed volume is appreciably 
diminished. The total reactor volume at 1 atm is 4060cm? 
(Table 5). By varying the pressure design to 4atm the 
total volume computed is 1960cm? for the same design 
conditions. The volume is significantly diminished (around 
52%). 

However, the kinetic parameters were obtained at atmo- 
spheric conditions; so, this behavior could not be exact. The 
pressure dependence of the kinetic rate expression used in the 
model needs further verification. In order to obtain more real- 
istic results it is necessary to determine kinetic data at different 
operation pressure values. 


4. Conclusions 


The results of this work show how model-based optimiza- 
tion at the design stage allows estimating efficiently the unit 
size of the WGS reactor, which is a critical component of 
ethanol processors, and optimizing the operating conditions. The 
methodology here applied is robust and can be used in a wide 
range of design cases. 

The WGS reactor is the largest and heaviest component 
because the WGS reaction is relatively slow compared to the 
other reactions involved in the whole reforming process, and is 
inhibited at high temperatures due to the thermodynamic equi- 
librium. 

The heterogeneous model allows computing the optimal reac- 
tor length and diameter, the optimal catalyst particle diameter, 
as well as the inlet reactor temperature of gaseous mixture for 
adiabatic operation mode. 

At the conditions here investigated, the results show that, it 
is necessary more than one adiabatic catalytic bed in order to 
achieve the specified conversion with the catalyst temperature 
imposed. However, if at the design stage, sizing is performed 
computing the reactive catalytic bed together with the insulating 
material it is possible to reach the design goal in a single unit. 
A suitable insulation material thickness may reduce the reactor 
volume since that heat loss to surroundings favors the reaction 
conversion. 

Although the results of the proposed optimization model 
depend on the catalyst characteristics and input and out- 
put conditions, the methodology and preliminary results 
obtained from this work are expected to be useful for process 
design, optimization and control of commercial fuel proces- 
sors for producing and purifying hydrogen for PEM fuel 
cells. 
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Appendix A. Estimation methods for gas mixture 
properties 


A.l. Viscosity 


The Bromley and Wilke modification of the theoretical 
Hirschfleder method was used to estimate the pure component 
viscosity (see Perry and Chilton [20]): 


33.3./ Mito; 
we (Losers 
Vc; 


0.261 
a .9Trp?? log; (1.9Tr;) 


The gas mixture viscosity was calculated from Wilke’s method 
(Reid et al. [19]): 


-1 


m=). Viki X yiii 
j 


i 


where gij; = pi! = (M\M;—')'?, as given by Herning and Zip- 
perer (Reid et al. [19]). 


A.2. Thermal conductivity 


The pure component thermal conductivities were calculated 
using the Eucken’s approximantion (Perry and Chilton [20]): 


—4 =f 2.48 
Aj = ui4.1890 x 10 Cp;2.394 x 107° + va 
i 
The following expression was used for computing the mixture 
thermal conductivity (Perry and Chilton [20]), 


1/3 
= Sr iyiki My! 
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A.3. Heat capacity 


The fluid mass heat capacity is calculated as 


where the pure component heat capacities (Cp;) were obtained 
from Reid et al. [19]. 


Appendix B. Insulation model 


In the thermally insulated reactor case, a heat loss towards 
the outside is considered modeling the energy transfer by con- 
duction through the insulation material, and by convection and 
radiation from the insulation surface to the environment. 

In Eq. (2), Q™ represents the heat flow per unitary reactor bed 
volume, and is related with the fluid temperature and the tube 
wall temperature by the heat transfer coefficient (hw) determined 
by the Leva correlation (see Froment and Bischoff [22]) 


mx 4 in 
Q = pho Tw) (B.1) 
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By applying the Fourier law to a cylindrical geometry and con- 
sidering that the thermal conductivity varies with temperature, 
the heat flow by conduction through the insulation material 
results in the following expression 

o™ 


k(T)dT 


i, 
(B.2) 
T. 


wo. 
DÈ In(Dins/ Di) 
Finally, the heat flow is related to the heat loss from insulation 


surface by convection (B.4) and radiation (Eq. (B.5)) (see Koenig 
[26]). 


ins 


AD ing 
O™ = — = (ev + qrad) (B.3) 
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1 
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qrad = F€ins((Tins + 273.16)* — (Tamb + 273.16)*) (B.5) 
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